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Abstract 

This paper is devoted to investigate f(R) gravity using Noether 
symmetry approach. For this purpose, we consider Friedmann Robertson- 
Walker (FRW) universe and spherically symmetric spacetimes. The 
Noether symmetry generators are evaluated for some specific choice 
of f(R) models in the presence of gauge term. Further, we calculate 
the corresponding conserved quantities in each case. Moreover, the 
importance and stability criteria of these models are discussed. 
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1 Introduction 



Astrophysical data from different sources such as cosmic microwave back- 
ground fluctuations [1J, Supernovae la (SNIa) [2] experiments, X-ray experi- 
ments [3] and large scale structure [I] indicate that our universe is currently 
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expanding with an accelerated rate. Higher dimensional theories [5] like M- 
theory or string theory may explain this accelerated expansion. Another 
explanation comes from modification of Einstein's theory with some inverse 
curvature terms which cause increase in gravity [6] . However, modified grav- 
ity with inverse curvature terms is known to be unstable and do not pass 
some solar system tests [7j. This discrepancy can be removed by including 
higher derivative terms. In particular, the viability jS] can be achieved with 
squared curvature terms. It is now thought that the current cosmic expan- 
sion can be justified if some suitable powers of curvature are added to the 
usual Einstein-Hilbert action [9]. Thus it would be interesting to investigate 
the universe in the context of modified or alternative theories of gravity. The 
f{R) theory of gravity, which involves a generic function of Ricci scalar in 
standard Einstein-Hilbert lagrangian, is an attractive choice. 

In recent years, many authors investigated f(R) gravity in different con- 
texts. Felice and Tsujikawa [10] gave a detailed review about f(R) theories 
of gravity. A similar work has been reported by Sotiriou and Faraoni [TT] . 
Hendi and Momeni [T2] explored black hole solutions in f(R) gravity. Moon 
and Myung [13J gave the stability analysis of the Schwarzschild black hole in 
this theory. Geodesic deviation equation in metric f(R) gravity is obtained 
by Guarnizo et al. [H]. Upadhye and Hu [15] investigated the existence 
of relativistic stars in f(R) gravity. The metric f(R) theories of gravity are 
generalized to five dimensional spacetimes by Huang et al. [TB] . They showed 
that expansion and contraction of the extra dimension prescribed a smooth 
transition from deceleration phase to acceleration phase. The stability con- 
ditions for f(R) models have been discussed by Starobinsky [17]. Multamaki 
and Vilja [18} IT9] explored spherically symmetric vacuum and non- vacuum 
solutions in f(R) theory of gravity. Shojai and Shojai [20] calculated exact 
spherically symmetric interior solutions in metric version of f(R) gravita- 
tional theory. Azadi et al. [21] investigated cylindrically symmetric vacuum 
solutions in this theory. Plane symmetric solutions are studied by Sharif and 
Shamir |22j . The same authors [231 121] investigated the solutions for Bianchi 
types / and V models for both vacuum and non-vacuum case. 

The field equations in f(R) gravity are fourth order partial differential 
equations (PDEs) when the function is assumed to have the terms like R 2 . 
However, the order could be higher if the terms like R 3 , R 4 etc. are included. 
On the other hand, Lie's theory gives a systematic and mathematical way 
to investigate the solutions of differential equations. The application of Lie 
group theory for the solution of nonlinear ordinary differential equation is one 
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of the most fascinating and significant area of research. It is mentioned here 
that from Lie's theory one can not only construct a class of exact solutions 
but can also find new solutions using different invariant transformations. It 
also gives most widely applicable technique to find the closed form solution 
of differential equations. Investigation of these solutions plays a vital role for 
the understanding of the physical aspects of these differential equations. 

Noether symmetry approach is an important aspect of Lie theory. This 
is the most elegant and systematic approach to compute conserved vectors, 
given by Noether in 1918. The conservation laws play a vital role in the 
study of physical phenomenon. The integrability for PDEs depends on num- 
ber of conservation laws. Another important aspect of conservation laws is 
that they are helpful in the numerical integration of PDEs. There are num- 
ber of methods developed for the construction of conservation laws such as 
Noether theorem [251 12S] for variational problems, partial Noether theorem 
for variational and non variational structures [27] and multiplier approach 
|28j . Computer packages for the construction of conserved quantities are 
reported by many authors. Some of them are: Wolf [21], Wolf et al. [30J, 
Goktas et al. [3 1J and Hereman et al. [321 1331 El]- The Maple code to 
compute conservation laws based on multipliers approach was introduced by 
Cheviakov [35J. 

Noether theorem states that any different iable symmetry of the action of 
a physical system has a corresponding conservation law. The main feature of 
this theorem is that it may provide information regarding the conservation 
laws in theory of relativity. Conservation laws of linear and angular momen- 
tum can be well explained by the translational and rotational symmetries 
using Noether theorem [36] . It has many applications in theoretical physics. 
In recent years, many authors have used this theorem in different cosmolog- 
ical contexts. Capozziello et al. [31] discussed f(R) gravity for spherically 
symmetric spacetime using Noether symmetry. Flat FRW universe has been 
discussed in Palatini f(R) gravity by Kucukakca and Camci [38]. Jamil et al. 
[39] investigated f{R) Tachyon model via Noether symmetry approach. Hus- 
sain et al. [30] found Noether symmetries for flat FRW model using gauge 
term in metric f{R) gravity. Energy distribution of Bardeen model is given 
by Sharif and Waheed [41] using approximate symmetry method. The same 
authors [42J re-scaled the energy in the stringy charged black hole solutions 
using approximate symmetries. In a recent paper, we [H] have found a new 
class of plane symmetric solutions in metric f{R) gravity using Lie point 
symmetries. 
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In this paper, we focus our attention to investigate the Noether symme- 
tries of FRW and spherically symmetric spacetimes in the context of metric 
f(R) gravity. The paper is organized as follows: In Section 2, we present 
some basics of f(R) theory of gravity. Sections 3 and 4 are used to cal- 
culate Noether symmetries of FRW and spherically symmetric spacetimes 
respectively. In the last section, we summarize the results. 



2 f(R) Gravity and Field Equations 

The action for four dimensional f(R) theory of gravity in gravitational units 
(87rG = 1) is given by [3j 



SfiR) = / V=g(f(R) + L m )d 4 x, (1) 



where L m is the matter Lagrangian and f(R) is a general function of the 
Ricci scalar. The standard Einstein-Hilbert action can be obtained by taking 
f(R) = R. Variation of this action with respect to the metric tensor yields 
the field equations 

f\R)R^ - \f(R)g»» - V,V v f\R) + g^Df(R) = kT^, (2) 

where prime denotes derivative with respect to R, n is a coupling constant in 
gravitational units, T™ is the standard standard matter energy-momentum 
tensor and 

□ = V% (3) 

with V M defined as the covariant derivative. The field equations can be 
expressed in an alternative form familiar with General Relativity (GR) field 
equations as 

1 

T 

where T™ u = T™/f'(R) and energy-momentum tensor for gravitational fluid 
is given by 



G ^ u R^iu n 9^uR + T^ v , (4) 



T ^ f'(R) 



\g»v(f{R) - Rf'(R)) + f\R)' aP {g a »gpu - g, v g aP )\ (5) 



It is clear from Eq.(T4]) that energy- momentum tensor for gravitational fluid 
T° contributes matter part from geometric origin. This approach seems 
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interesting as it may provide all the matter components which are required 
to investigate the dark part of our universe. Thus it is expected that f(R) 
theory of gravity may give fruitful results to understand the phenomenon of 
expansion of universe. 



3 Noether FRW Symmetries 

In this section, we shall find the Noether symmetries of FRW spacetime. The 
FRW metric is given by 

7 2 7 2 2 df 

as — at — a 



where a is function of cosmic time t and known as scale factor of universe 
and dVL 2 = r 2 {dQ 2 + sin 2 9d<fi 2 ) . The curvature parameter k is 0, 1 or — 1, 
which represents flat, open or closed universe respectively. The corresponding 
Lagrangian is given by jH] 

L = Qad 2 f + Qa 2 Raf" + a 3 (f - Rf) - Qkaf + a 3 P. (7) 

Here dot denotes derivative with respect to t and the fluid pressure P is given 
by 

P = mwa- 3(1+u) , (8) 

where m is an arbitrary real constant while u is the equation of state param- 
eter. Noether symmetry generator of Eq.([7]) is given by 

d d d 

X = r(t, a,R)-^ + 1>(t, a > R )^+ 0(*> a ' R )qr ■ ( 9 ) 

We can find Noether symmetries by using the following equation: 

X [1] L + (Dt)L = DB(t,a,R). (10) 

where X^ is the first prolongation [15] given by 

XW = X + a,R)— + a, R)-r, (11) 



where 



■dip dib dr -dib dr -dr , , 

• d(p d(j) -d(j) -dr - 2 dr -dr . . 

^=m +a da- + R dR- R m- R dR- aR d-a {U) 
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B is called the gauge function with D defined as 

d d ■ d 
D =M +a 3-a +R a-R- < 14 > 

The first integral is also known as conserved quantity associated with X and 
is defined as 

I = B - tL - (V - rd)^ - (<j> - tR)^. (15) 



i.e. 



Using (J7D in Eq. fllQI) . we obtain an over determined system of linear PDEs, 

Ta = 0, (16) 

TR = 0, (17) 

/>« = 0, (18) 

6a 2 /'ty t = B R (19) 

12a/V t + 6a 2 /'Vt = £«,(20) 

iPf' + (j)af" + 2af'<ip a -af'T t + a 2 f" ( j) a = 0, (21) 

2afil> + a 2 f M <j> + a 2 f' , (il> a + <j> R -T t ) + 2af , il> R = 0, (22) 
a 2 (3ip + ar t )(f - Rf) - a^Rfcj) + T t (-6kaf + wmcT 3w ) - 

tfj(6kf + 3u 2 ma- 3uJ - 1 ) = B t .(23) 



We solve these equations using different assumptions. From Eq.f ll8l) . we have 
f" = with ip R 7^ 0. This further gives /' = by using Eq. (j220 . Solving 
determining equations with these conditions, we obtain trivial solution, i.e. 

r = 0, ^ = 0, = 0, (24) 

with zero gauge term. Hence we investigate the solution of determining equa- 
tions by taking ip R = with f" ^ 0. The solutions are discussed mainly for 
two different cases, namely flat vacuum universe and non-flat non-vacuum 
universe. 

Case 1: Flat Vacuum Universe (k = = tu) 

Here we solve the determining equations for flat vacuum case. It is mentioned 
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here that this case has already been discussed by Jamil et. al [40 J but ob- 
tained results suggest that the gauge term is zero. However, we have explored 
a more general solution of the determining equations for f(R) = foR 3 ^ 2 which 
gives a non-zero gauge term. The solution in this case is given by 



Cit + c 2 , 

2c\a 2 + 3c 3 t + 3c4 
3a 

Ci<2 2 + C 3 t + C 4 
— Zrt . 



T 

J, 



B = 9c 3 aVR + c 5 . 
The Noether symmetry generators turn out to be 

x 2 
x 3 
x 4 



t^ + ~a^-2R±, 



2 d 
-a— 

3 da 



0_ 

dt 
d_ 



(25) 
(26) 

(27) 
(28) 

(29) 
(30) 
(31) 
(32) 



These generators form a four dimensional algebra with the following commu- 
tator table: 



Table 1: Commutator Table 




x x 


x 2 


x 3 


x. 


x 1 





x 2 


3 


4 a: 4 

3 


x 2 


-x 2 





-x 4 





X 3 




x 4 








x, 


3 












where Lie bracket [Xi, Xj] is defined by the following unique relation 



[X h Xj] 



X.; ( Xi ) — Xa ( X 



where i,j = 1, 2, 3,4. Moreover, the first integrals in this case are 



h 
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9ata z R^ + -aHB? + -a'HaRR^ - 3a 2 aR^ - 3a 4 RR^ 



i i 9 • -i 

h = 9ai?5 _ 9taRz - -taRR~ , 

h = -9d#5 - - a RR^-. 

2 

Case 2: Non-Flat Non- Vacuum Universe (k ^ 0, w^O) 
For this case, the determining equations yield a solution 

r = ci, £ = 02, (33) 

^ = 0, = 0. (34) 

Here the gauge term turns out to be constant which can be taken zero. It 
is mentioned here that this solution is for an arbitrary f(R). The Noether 
symmetry generator turns out to be 

X = §. (35) 

The corresponding first integral becomes 

I = Qaa 2 f + Qa 2 dRf" - a 3 (f - Rf) + Qkaf - mujaT^ . (36) 

4 Noether Symmetries of Spherically Sym- 
metric Spacetime 

In this section, we shall find the Noether symmetries of static spherically 
symmetric spacetime 



1 «.„a 9 



ds 2 = Adt 1 



2 



— + r 2 dn 2 

A 



(37) 



where dfl 2 = r 2 (d9 2 + sin 2 9d(f>) and A is the function of r. The corresponding 
Lagrangian is given by 

L = ^-*n + 2f(l-A-r%) + r^)(£), (38) 
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here prime denotes derivative with respect to R. 



Corresponding Noether symmetry generator is given by 

X = r(r, + V(r, R, A)-^ + <f>(r, R, A)-^. (39) 

The Noether symmetries can be computed by using the following equation 

X [1] L + (Dr)L = DB(t,a,R) } (40) 
where is the first prolongation [45j given by 

XW = X + 4(r, R,A)-r + 0(r, R, A)—r, (41) 

oR oA 

in which ' represents derivative with respect to r and 

' dip A dib >8t >dib > 9 dr > > dr , 

* = w +R m- R * +A iix-*m-™M> (42) 

V dr dR dA dr OA OR K J 

In Eq.( l40l . B is called the gauge function with D defined as 

d d ■ d 

Substituting ( 138]) in Eq. (140|) and after some manipulations, we get an over 
determined system of linear PDEs, i.e. 







= o, 


(45) 






= o, 


(46) 




^4 


= o, 


(47) 




<pR 


= o, 


(48) 






= Br, 


(49) 


2r( 


rf" + fr r ) + rW + 0a/" + M" - r r f) 


= o, 


(50) 




-2rf'-2r^f" + f ( t> A )+r 2 f"^ r 


= Ba, 


(51) 


2rr(/ 


- Rf) - 2 ftp - r 2 Rf"iP + 2/> - 2/>A - 








2rf<p r + r 2 r r / - r 2 r r i?/' + 2r r f - 2r r Af 


= B r . 


(52) 
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When f(R) is arbitrary, we obtain trivial results, i.e. 



r = 0, ^ = 0, = 0, (53) 

However we use a well known form of f{R), i.e. 

f(R) = f R n , n ± 0, 1. (54) 

This function has been widely used in different cosmological context. The 
determining equations, using Eq. (l5"4j) yields 

R c 

T = Cl r, ^ = -3ci— (p= — (2An-3A-2n + 3), (55) 
n n 

where the gauge term is zero in this case. Thus the Noether symmetry 
generator in this case is given by 

X = r£ - 3- A + \ 2 An -ZA-2n + 3) A (56) 
ar n ait n <M 

and the conserved quantity turns out to be 

J = r(n-l) [r 2 /? n +6(A-l)i?"- 1 +5rii? n - 1 -(2n-3)(A-l)r J R J R n - 2 +r 2 ii? J R n - 2 

(57) 



5 Summary and Conclusion 

The main objective of this paper is to investigate the Noether symmetries in 
metric f(R) gravity. FRW and spherically symmetric spacetimes are consid- 
ered for this purpose. We present a general solution of determining equations 
for FRW universe with gauge term. In fact, we get four Noether symmetry 
generators. A non-zero gauge term is obtained which depend on Ricci scalar 
R and scale factor a. It would be worthwhile to mention here that solution 
already obtained by Jamil et al. jlQ] is a subcase for the flat universe with 
zero gauge term. Moreover, in palatini f(R) gravity, a non-zero time de- 
pendent gauge term is obtained [SB]. For non-flat universe, we obtain one 
symmetry generator which is translation of time coordinate. The first in- 
tegrals are obtained in each case. The interesting feature of cosmological 
model for f(R) = f R 3 / 2 is that it gives a negative deceleration parameter 
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which is consistent with the recent experimental results to justify the accel- 
erated expansion of universe. It has been shown [H] that for a = dot 2 and 
f(R) = foR 3 ^ 2 , the deceleration parameter is — \ for the flat universe. 

The spherically symmetric spacetimes yields a set of eight linear PDEs. These 
equations are solved for two cases of f{R)~- First case involves an arbitrary 
function of Ricci scalar which gives a trivial solution. However, we obtain a 
non-trivial solution in second case when f(R) = foR n - The corresponding 
conserved quantity is also obtained in this case. This cosmological model has 
been used extensively in the recent literature. In particular, the well known 
f(R) model with inverse curvature term, corresponding to n = —1, predicts 
late time accelerated expansion of the universe [17] . 

Moreover, the stability conditions for f(R) models are f'(R) > and f"(R) > 
[IT]. It would be worthwhile to mention here that the model f(R) = foR 3 ^ 2 
satisfy these conditions for f Q > and R > 0. These conditions are also sat- 
isfied by the model f(R) = f R n when f > 0, n - 1 > and R > 0. 
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